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AMnct— An  unstable  pitch  control  systeirTllf  a  terminal  -homing  missile  was  formerly  stabilized  using  a 
high  order  stabilization  filter  that  was  realized  using  active  elements.  A  new  dominant-data  matching 
method  is  presented  to  redesign  the  high-order  stabilization  filter  for  obtaining  reduced-order  filters.  As  a 
result,  the  implementation  cost  is  reduced  and  the  reliability  increased.  An  algebraic  method  is  also  applied 
to  improve  the  performance  of  the  redesigned  pitch  control  system.  In  addition,  the  proposed  dominant- 
data  matching  method  can  be  applied  to  determine  a  reduced-order  model  of  a  high-order  system.  Unlike 
most  existing  model  reduction  methods,  the  reduced-order  model  has  the  exact  assigned  frequency-domain 
specifications  of  the  original  system.  Computer-aided  design  methods  can  also  be  applied  to  design  general 
control  systems. 


I.  INTRODUCTION 

The  pitch  control  system  of  an  unstable  terminal  homing  missile  [1]  was  formerly  stabilized 
using  a  fourth-order  series  compensator.  The  compensator  had  two  pairs  of  complex  poles  that 
were  realized  using  active  elements.  The  objective  of  this  paper  is  to  develop  a  computer-aided 
method  for  redesigning  the  compensator  such  that  the  implementation  cost  of  the  compensator 
can  be  reduced  and  the  performance  of  the  redesigned  pitch  control  system  improved.  The 
block  diagram  of  the  existing  stabilized  system  is  shown  in  Fig.  I,  and  its  over-all  transfer 
function  is 


■r-/  ,  T,(s)Gp(s)  _  G,(s) 

I  +  T,(s)G<ts)H,(s)  l  +  c,(j) 

_  bp+  bis  +  •  •  •  +  b9s*  +  hl0s"  a  N{s) 
a0+a,s+  •  •  ■  +aioJlo+aiiJn  ~  D(s) 


where 


a*  =8.802158509x  10" 
a ,  =  2.419047424  x  10" 
a2  =  2.911920560x  10" 
a,  =2.420405431  x  10" 
a4  =  6.667397031  x  10" 
a,  =  9.749923212  x  10" 
=  9.360329977  x  10*1 
u,  =6.231675318x  10" 
a ,■  =  2.9769506%  x  10* 
a,  =  9.316239040  xIO5 
a,«=  1 .923554000  x  103 
«„  =  I. 


b0  =8.80215809x  10" 
h,  =4.610004670x10" 
b2  =  2.926344345  x  10" 
by  =5.017212044x  10" 
bt  =  2.563396371  x  10" 
by  =1.494523312x  10" 
bo  =0. 


bn  =0. 
b»  =0. 
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Fig.  1.  The  block  diagram  ot  the  existing  control  system. 


and 


T,(i)  =  The  transfer  function  of  the  existing  series  stabilization  filler 


*  N,(s) 

D,U) 


(1c) 


Hf(i )  =  The  transfer  function  of  the  gyro  =  I  (Note  thal  a  rate  gyro  is  not  available)  (Id) 

Go(j)  =  The  transfer  function  of  the  actuator  and  air  frame  dynamics  of  the  missile 
system— The  plane  transfer  function 

*  The  open-loop  transfer  function  of  the  original  pitch  control  system  with  T,{s)  =  I 
and  H,(s)  =  I 

324332.3I6Q  +  O.I933Xi  +  63Xs  +  1500) 

*  id  -  2.920(1  +  3.173X1  +  87.9*  /95.5Xi  +  ll2.3Hi  +  1383)  '  ' 

G,(s)  =  The  open-loop  transfer  function  of  the  existing  stabilized  system 

-  T,(i)G«(i)W,(i).  (If) 

For  ease  of  presentation,  we  define  the  dominant  data 


(a)  The  real  and  imaginary  parts  of  the  transfer  function  when  <v  =  0. 

(b)  The  gain  margin. 

(c)  The  phase-crossover  frequency. 

(d)  The  phase  margin. 

(e)  The  gain  crossover  frequency. 

Nyquist  plots  of  G,(i)  and  Gg(i)  are  shown  in  Fig.  2.  The  dominant  frequency-response  data  of 
G,(i)  are  defined  as  follows 


(I)  The  real  (ft,)  and  imaginary  (fm)  parts  of  G,(.r)  at  i  =  jtai  =  jO  are 


R,GAj 0)  =  -  2. 103817  and  /„G,</0)  =  * 
(2)  The  gain  margin  Gm  of  this  system  is 


i 


G,(/o»„)|  \R,G.{h>n) 


I 


1.5 


(2a) 


(2b) 


where  the  phase-crossover  frequency  is 


L 


G,(fu>„)  =  -  180“  or  u„  -  l.9rad./sec. 


(2c) 


in 
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The  equivalent  real  and  imaginary  parts  of  G,(jo>„)  at  u>„  =  1.9  are 


R,G,(m.)  =  - 1.507944  (2d) 

LG,(h>„)  =  -  0.006490205.  (2e) 

(3)  The  phase  margin  <t>m  of  this  system  is 

4>„  =  180°  -t-/ G,(j(o„)  a  5.7787°  (2f) 

where  the  gain-crossover  frequency  <oK  is 

|G,(/«w)|  =  1  or  av  s  3.2  rad/sec.  (2g) 

The  equivalent  real  and  imaginary  parts  of  G,(ja>K)  at  wK  =  3.2  are 

U,G,(/w„.)  =  -0.9939143  (2h) 

lmG,{ja>K)  =  -  0.09547478.  (2i) 


The  frequency-response  data  at  <o  =  0  in  (2a)  indirectly  indicates  the  steady-state  value  of  the 
unit-step  response  of  T,(s).  The  data  at  to  =  o>„  and  <o  =  aw  in  (2)  represent  two  control 
specifications  (2]:  gain  margin  and  phase  margin  that  characterize  the  relative  stability  and  the 
transient  response  of  the  existing  stabilized  system.  The  dominant  frequency  response  data  of 
the  plant  G^s)  and  others  are  as  follows 

(I)  The  real  and  imaginary  parts  of  GoO'at)  at  o>  =  0  are 


R,G„(/0)  =  - 1.304841  and  ImG<jj0)  =  « 

(2)  The  phase  margin  fom  of  the  plant  is 

<frv,=  180°  +  /Gq(mJ=-  5.58° 

where  the  gain-crossover  frequency  an*-  is 

|Go(/woc)|  =  I  or  a**,  a  1.6  rad/sec. 

Other  frequency-response  data  at  a>„  and  aw  are 

R,Go(ja>„)=  -0.93707661  . 
fmGo(/a>,»)  =  0.06716120  }  for  ~  1  9 


R,G<Aja>")  =  -0.6181657]  . 

(4)  Z.GoOw.)  =  0.01949691  }  forw"-3'2' 


(3a) 


(3b) 


(3c) 


(3d) 


(3e) 


Comparing  the  dominant  data  of  G,(s)  in  (2)  and  G^s)  in  (3d)  and  (3e)  yields  the  dominant 
data  required  of  the  stabilization  filter  T,(s)  as 

(1)  RrT,(jat)  =  1.6  and  lmT,(ja>)  =  0  at  w  =  0  (4a) 

(2)  R,T,(ftt„)  =  1.600492  and  /.!,(/»„)  =  0.1216316  at  =  1.9  (4b) 

or 

|T,(/ww)|  =  1.605107127  and  /r,(K,)  =  4.345918I9T  at  =  19  (4c) 

(3)  R,T, (/aw)  =  1.601402  and  lmT,{j»«)  *  0.2049554  at  «w  -  3.2  (4d) 
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|T,(jtaw)|  =  1.614464333  and  /r<(/n»,)  =  7.293349493°  af  =  3.2.  (4e) 

Fiona  (le)  and  (10  we  observe  that  G^s)  and  G,(s)  are  non-minimum  phase  functions.  From 
the  Nyquist  plots  in  Fig.  2  and  the  Nyquist  stability  criterion,  we  can  conclude  that:  the  original 
missile  system  (without  the  stabilization  filter)  is  unstable,  the  existing  stabilized  system  is 
asymptotically  stable  and  its  time  response  is  oscillatory  due  to  the  small  positive  phase  margin 
in  (20-  The  purpose  of  this  paper  is  to  develop  computer-aided  design  methods  for  redesigning 
the  stabilization  filter  to  reduce  the  implementation  cost  and  improve  the  flight  control 
performance  of  the  missile  system. 

Two  computer-aided  methods  are  developed  in  this  paper  and  subsequently  used  to  redesign 
the  pitch  control  system.  In  Section  II  a  dominant-data  matching  method  for  modeling  a 
transfer  function  (called  a  standard  transfer  function  Tr(s))  that  matches  the  assigned 
specifications  shown  in  (2)  is  developed.  The  obtained  standard  transfer  function  T,{s )  is  a 
reduced-order  model  of  the  existing  stabilized  system  T,(j)  in  (la).  The  time  and  frequency- 
response  curves  of  T,{s)  and  T,(s)  will  be  compared  to  verify  that  the  data  in  (2)  are  dominant. 
The  dominant-data  matching  method  is  then  applied  to  obtain  the  reduced-order  model  of  the 
existing  stabilization  filter  T,(s).  Also,  the  method  is  used  to  fit  a  low-order  model  that  satisfies 
the  specifications  shown  in  (4).  Thus,  two  low-order  stabilization  filters  are  obtained.  In  Section 
III,  we  apply  the  dominant-data  matching  method  and  the  algebraic  method  due  to  Shieh[3]  and 
Chen[4)  to  redesign  the  pitch  control  system.  In  order  to  simplify  the  design  process,  the 
dominant-data  matching  method  is  first  applied  to  obtain  an  unstable  reduced-order  model  of 
the  original  unstable  high-order  system  G^s).  Then,  the  algebraic  method  is  applied  to  redesign 
the  pitch  control  system  that  has  a  series  filter  in  the  feed  forward  loop  and  a  parallel  filter  in 
the  feedback  loop.  Thus,  the  advantages  of  a  compensator  in  the  feedback  structure  can  be 
fully  used. 

II.  THE  DOMINANT-DATA  MATCHING  METHOD 

The  design  goals  and/or  the  nature  of  the  transient  response  of  a  control  system  are  often 
characterized  by  a  set  of  control  specifications  [2].  These  specifications  are  commonly  classified 
as:  (I)  time-domain  specifications,  e.g.  rise  time  and  overshoot;  (2)  frequency-domain 
specifications,  e.g.  gain  and  phase  margin;  (3)  complex-domain  specifications.  e.g.  damping  ratio 
and  underdamped  natural  angular  frequency.  Rules  of  thumb  that  represent  the  relationships 
among  the  above  three  control  specifications  have  been  proposed  by  Axelby  [5]  and  Seshadri  ft 
of.  [6].  Given  the  above  guidelines,  it  is  obvious  that  the  gain  margin,  phase  margin,  phase- 
crossover  frequency  and  gain-crossover  frequency  are  the  most  important  specifications.  These 
data  are  called  the  dominant  frequency-response  data  and  shown  in  (2b)  through  (2i).  Other 
important  frequency-response  data  is  the  steady-state  value  of  a  closed-loop  system  that  is 
indirectly  represented  by  the  value  of  G,(j» )  at  w  *  0  in  (2a).  These  dominant  data  in  (2)  may  be 
considered  as  the  design  goal.  In  order  to  verify  that  the  data  in  (2)  are  dominant  ones  we  need 
to  construct  a  transfer  function  T,(s)  that  is  the  reduced-order  model  of  T,(j)  and  has  the  exact 
assigned  specifications  shown  in  (2).  The  time  response  curve  and  the  corresponding  time- 
domain  specifications  of  this  practical  system  T,(s )  are  difficult  to  obtain  because  T,(s)  is  a 
high-order  transfer  function  with  large  coefficients.  Furthermore,  it  is  a  stiff  function.  The  latter 
can  be  verified  from  its  small  coefficient  aM  (the  sum  of  all  poles  of  the  system)  and  its  large 
constant  u«  (the  product  of  all  poles  of  the  system)  in  (la).  As  a  result,  many  numerical 
integration  methods  (for  example,  the  Runge-Kutta  method  (7|)  require  large  data  size  and  high 
precision  calculation  for  time  response  determination.  However,  the  frequency  response  curve 
and  the  corresponding  frequency-domain  specifications  of  this  system  can  be  easily  determined 
by  a  digital  computer.  Thus,  a  frequency-domain  approach  or  a  dominant-data  matching  method 
is  proposed  to  determine  the  reduced-order  model  T,(s)  and  to  redesign  the  pitch  control 
system.  There  exists  several  frequency-domain  methods  for  model  reductkms(S-l01.  However, 
the  reduced-order  models  obtained  from  the  proposed  method  gives  the  exact  assigned 
frequency-domain  specifications  in  contrast  to  other  techniques  which  do  not. 

Let  tile  desired  reduced-order  model  of  T,(s)  or  the  standard  model  that  represents  the 
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design  goal  in  (2)  have  the  following  form 


Ti  >  ao  +  biS  +  bjS2  GM 

r's  ~  ao+ a, s  +  a:s*+ s’  I  +  Gr(s) 


where  the  open-loop  transfer  function  G,(r)  is 


.  ao  +  bjS  +  biS 

''  ~  rKai-hil  +  lfli-hjIs  +  s4!' 


(5a) 


(5b) 


The  unknown  constants  u,  and  6,  are  to  be  determined  from  the  conditions  in  (2).  Following  the 
basic  definitions  and  knowing  the  required  values  from  (2)  yields  a  set  of  nonlinear  equations 
/((a#,  a,.  a2,  6,,  b2)  =  0  for  i  =  1, 2 . 5  as  follows 


(I)  The  requirement  of  (2a),  or  R,G,(/0)  =  -  2.1,  gives 


/,(flo.  a,,  a2.  6,,  b2 )  =  a,b,  -  6,*  -  a„a:  +  a<,6j  +  2.1(at  -  6,)*  =  0  (fca) 

(2)  The  requirement  of  (2b),  or  R,G,(/w„)  =  -1.5  at  =  1.9  gives 

/}(do.  d|,  flj,  6,,  6j) =  (flj  —  hjXu#  -  3.6162) —  6, (a,  -  hi  —  3.61) 

-  l.5(3.6l(aj-  hj)}  +  (a,  -  6,  -  3.61)*]  =  0  (6b) 

(3)  The  requirement  in  (2c),  or  / Gr(jai„)  =  -  180°  at  u>„  =  1.9,  gives 


/](<io,  ii|,  di,  b  i,  6j)  —  3.61h|(fl2  -  6>)  +  (o«  —  3.6I6jX<Ji  —  6,  —  3.61)  =  0  (6c) 


(4)  The  requirement  of  (If),  or  4*.  =  5.7®  at  =  3.2,  yields 


/<(uo.  <ii.  Uj.  6,,  6j) =  I0.24h|(aj  -  hj)  +  (a# -  I0.246jXa,  -  6|  -  10.24) 
-0.3l040224{(fl3-  hjXflo-  I0246:) 

-Mu, -6,-10.24)1  =  0  (6d) 

(5)  The  requirement  of  (2g),  or  |G,(jw„.)|  *  1  at  uK  =  3.2,  gives 


/da#,«„a2,h„hJ)  =  (ai,-  10.246,)*  +  10.246,*-  t(»4.8576(u,  -  6,)* 

-  I0.24(a,  -  6,  -  10.24)*  =  0.  (6e) 


The  above  set  of  high-order,  nonlinear,  algebraic  simultaneous  eqns  (6)  are  difficult  to  solve. 
The  Newton-Raphson  method [II]  that  is  available  as  a  library  computer  program  package 
(called  the  Z  system!  12]  in  many  digital  computers  can  be  used  to  solve  these  nonlinear 
equations.  However,  it  is  well  known  that  the  Newton-Raphson  method  will  only  converge  to  a 
desired  solution  for  a  small  range  of  starting  values  or  initial  estimates.  In  order  to  improve  the 
convergent  speed,  the  following  methods  are  suggested  for  good  initial  estimates. 

(I)  Initial  estimate  by  using  the  model  reduction  method  due  to  Shieh  and  Chen[3,8]. 

Shieh[3)  and  Chen  [8]  have  proposed  a  continued  fraction  method  for  model  reduction.  The 
method  is  as  follows.  The  N(s)  and  D(s)  in  (la)  are  arrayed  into  ascending  order  and  expanded 
into  the  continued  fraction  of  the  second  Cauer  form  by  performing  repeated  long  divisions.  i.e. 


T.(s) 


N(s)  ^  6p  +  6,s  +  ■  •  •  +  6,0*“ 
D(i)  *  <4  +  «,5  +  •  •  •  +  a,,j" 


*■  +  ■ 


hj+- 


*i  +  - 


6«  +  - 


(7a) 


*  - 
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where  *,-1.  *2  = -0.401749,  *,=  -0.475321,  *4  =  25.1998,  *,  =  -0.0322195,  At  = -24.1061, 
*7=  •  ■  •,  •  •  •,  and  *22  =  ‘  • 1  for  the  example  problem.  The  reduced-order  can  be  obtained  by 
retention  of  the  first  several  dominant  quotients  (A,)  for  /  =  1, 2, ...  as  follows 


fit 

*1*2  +  5 


1 


*1  + 


5 

*2  +  sl  hy 


hihj  +  5 

*1*2*3  +  (*|  +  *3)5 


*2*3*4  +  (*2  +  *4)5 _ 

*,*2*3*4 +  (*1*2  +  *1*4+  *3*4)5  +SJ 

_ *2*3*4*,  +  (*2*3  *2*5  ^  *4*5)5  +  5* _ 

*  I  *2*3*4*,  +  (*1*2*3  +  *1*2*5  +  *1*4*3  +  *3*4*,)5  +  (*|  +  *3  +  *,)s2 

*2*3*4*, *6  +  (*2*3*4  ^  *2*3*6  +  *2*3*6  +  *4*5*615 

*  |  *2*3*4*,  *6  +  (*,*2*3*4  +  *,*2*3*6  +  *1*2*, *6  +  A,  *4*5*6  +  *3*4*, *6)5 

_ +  (*2  +  *4  +  hf,)s* _ _ _ 

+  (*,*2  +  *,*4  +  *,*6  +  *3*4  +  *3*6  +  A,  *6)5*  +  5* 


(7b) 


(7c) 


(7d) 

(7e) 


(7f) 


where  (7f)  is  the  3rd  order  approximate  model  of  the  original  11th  order  system,  or  in  our 
problem 


3.7376  +  19.46925  +  0.6920s* 

,{s'  3.7376+ 10.1661s  +  0.9488s*  +  s’  '  ' 

Using  the  coefficients  in  (8a)  as  initial  estimates:  aj  =  3.7376,  at  =  10.1661,  of  =  0.9488, 
*t  =  19.4692  and  *f  =  0.6920  and  applying  the  Newton-Raphson  method  [12]  to  solve  the 
nonlinear  equations  in  (6)  yields  the  desired  solution:  a«  =  6.37807,  a,  =  10.46222,  a2  =  1 .259008, 
A,  =  20.55661  and  *2  =  0.243466  at  the  8th  iteration  with  the  error  tolerance  of  10-6.  The  desired 
reduced-order  or  the  standard  model  is 

r  6.37807  +  20.55661  s  + 0.243466s*  0,(5) 

f,w  6.37807  +  10.462225  +  1 .259008^  +  ?  1  +  Gr(s) 


where 

G,(s )  =  The  open-loop  transfer  function  of  the  standard  model 

6.37807  +  20.55661 5  +  0.243466s2 
s<  -  10.09439+  1.015542s  +  s1)  ’ 

The  Nyqtnst  plot  of  G,(s)  is  shown  in  Fig.  2,  the  unit-step  responses  of  T,(s)  in  (la)  and  Tr{s)  in 
(8b)  are  shown  in  Fig.  3.  The  approximate  results  are  satisfactory.  Thus,  we  verify  that  the  data 
in  (2)  are  dominant.  It  is  obvious  that  although  the  T*,(s)  in  (8a)  may  have  a  good  overall 
approximation  of  T,(s),  only  T,(s)  in  (8b)  has  the  exact  assigned  frequency-domain 
specifications  as  required  in  (2)  that  is  essential  in  the  design  of  a  control  system  in  the 
frequency  domain.  It  may  also  be  noticed  that  the  T*r(s)  obtained  by  the  continued  fraction 
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Fig.  3.  Time  response  of  various  models. 


method  (7)  may  be  unstable  even  if  the  original  system  T,(s)  is  stable.  The  following  mixed 
method  is  suggested  for  obtaining  a  stable  reduced-order  model. 

(2)  Initial  guess  by  using  the  mixed  method. 

The  mixed  method  has  both  advantages  of  the  continued  fraction  method  [3, 8],  the 
dominant-pole  method!  13]  or  the  equivalent  dominant-pole  method [9],  It  can  be  applied  to 
determine  a  stable  reduced-order  model  from  which  a  good  initial  estimate  can  be  determined. 
The  mixed  method  is  as  follows.  The  relationship  between  the  quotients  A,  and  the  coefficients 
a,  and  b,  in  (7a)  can  be  expressed  in  the  following  matrix  equation(3,4]: 

l*]  =  ltfM  (9) 


where 


{a]r  =  {flo.<ii,<i: . a„-i] 

[fr]r  =  ffro,h„fr2 . h„-,] 

[H]  -Wttf,] 


and 


[f/2}  = 


A,  0  0 
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0  0  1 
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T  in  (9)  designates  transpose.  The  desired  reduced-order  model  may  be 


T,(s)  = 


d0  +  d\S  +  ■  •  •  + 

Zo  +  eis  +  •  •  •  +  *,-!$'''  + 


=  !. 


(10) 


The  coefficients  e,  in  (10)  can  be  determined  from  the  coefficients  of  the  polynomial  that  is  the 
product  of  the  dominant  poles  of  T,(s)  in  (7a).  Replacing  a,  in  (9)  by  e,  in  (10)  and  solving  the 
matrix  equation  in  (9)  yields  the  desired  coefficients  d,  in  (10).  The  obtained  Tr(s)  has  the 
dominant  poles  and  quotients  of  T,(s)  and  is  always  stable.  When  the  roots  of  D(s)  in  (7a)  are 
not  available,  the  approximate  equivalent  dominant  poles  and  the  resulting  coefficients  e,  can  be 
determined  from  the  Routh  table  as  suggested  by  Hutton  and  Fried!and[9].  Because  the  method 
uses  the  dominant  quotients  of  the  original  system  and  the  equivalent  dominant  poles  from  the 
Routh  table,  this  method  may  be  conveniently  called  the  mixed  method  of  the  continued 
fraction  approximation [8]  and  the  Routh  approximation [9],  Using  the  mixed  method  we  can 
obtain  another  stable  reduced-order  model.  Thus,  a  good  initial  estimate  can  be  determined. 

Since  the  transfer  function  of  the  existing  stabilization  filter  Ts(s)  is  available,  we  will  use 
the  proposed  dominant-data  matching  method  for  determining  the  reduced-order  of  T,(s).  The 
T,(s)  in  (lc)  can  be  considered  the  closed-loop  transfer  function  of  a  control  system 

r<  i  -Ns(*)  C.(s)  460800s2  +  69120000s +  144  x|Q7 

D,(s)  ~  1  +  G,(s)  s' + 250s5 + 76900s2 + 72 x  K>5s  +  9x  10"  Ula) 


where  the  open-loop  transfer  function  G,(s)  is 

r,  460800s2  +  69120000s  4- 144  x  IQ7 

0,'S'  si  +  250s5  -  383900s2  -  61920000s”-  J.4  x  10* 

The  dominant  frequency-response  data  of  this  system  are 


(lib) 


(1)  G,0'0)  =*  - 2.667 

(2)  R,G,(M,)  =  -  1.032833 
LGAju>")  =  0.002017351 

where  a>,„  =  the  phase-crossover  frequency  =  140  rad/sec. 

(3)  RflM+c)  =  ~  1002941 
ImGAk>,A  -  -0.03668759 

where  a>„  =  the  gain-crossover  frequency  =  200  rad/sec. 

The  reduced-order  stabilization  filter  T,,(s)  is  assumed  to  be 

t  t  ^  ^  b\5  _  Gs i(s) 

flo  +  fl|S  +  S^  1  4- G,|(s) 

where 


G.i($)  = 


_ b0  +  b,s _ 

(a0-b(i)  +  (ai-b,)s  +  s2 


(12a) 

(12b) 

(12c) 


(1 2d) 
(12e) 


(13a) 


(13b) 


and  a,  and  b,  are  unknown  constants  to  be  determined. 

Using  the  coefficients  at  and  b{  of  G,,(s)  and  following  the  basic  definitions  of  the  data 
shown  in  (12a),  (12b),  (12c)  and  (12d)  results  in  the  following  nonlinear  equations,  respectively 


(1)  G2b) 

a,.  bt)  *  l.6ad-0.6ao  -  19600)  +  I9600*,(a,  -  bt)  +  1.032833[(0.6ao 

+  19600)2  +  19600(O|  -  b|)2]  =  0  (14a) 
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(2)  (12c) 

Ma0,  a„  b,)  =  1406|(  -  0.6ao  -  19600)  -  224ao(a,  -  b,)  -  0.00201735 1 

[(0.6a0  +  19600)2  +  19600(0, -*,)2]  =  0  ( 14b) 

(3)  (12d) 

Mac.  a,,  b,)  =  1.6a0(  -0.6ao- 40000)+ 40000b, (a,  -  b,)+  1.002941 

l(0.6a,,  +  40000)J  +  40000(0,  - 1»,)2)  =  0  (14c) 

where  b0  =  I.6a0  as  obtained  from  (12a). 

The  initial  estimates  can  be  obtained  from  the  reduced-order  model  of  Tjs)  in  (1c)  using  the 
mixed  method  of  the  continued  fraction  approximation  and  the  Routh  approximation.  The 
reduced-order  model  is 

1281.40525s +  29937.62994 

r,{s>  r2  + 52.4375s  + 1871 1.01871’  '  ' 

Using  the  coefficients  in  T? ,(s)  as  initial  estimates  and  applying  the  Newton-Raphson 
method  [12]  we  have  the  desired  a,  and  b,  in  (13)  for  the  7th  iteration  with  an  error  tolerance  of 
10“*.  The  desired  low-order  stabilization  filter  is  therefore 

r  957.260014s  +  33467.93525 

s2  +  29.981293 s  +  20917.459536'  (  ' 

The  unit-step  response  of  the  existing  stabilized  pitch  control  system  in  (la)  and  the  redesigned 
pitch  control  system  using  T,,(s)  in  (16)  and  Go(s)  in  (le)  are  shown  in  Fig.  3.  The  result  is 
considered  to  be  satisfactory. 

An  alternate  approach  is  proposed  for  redesigning  the  stabilization  filter  as  follows.  Because 
the  function  of  a  stabilization  filter  is  to  convert  the  dominant  data  at  to  =  0,  to„  =  1.9  and 
b>K  =  3.2  of  the  original  unstable  system  G^s)  in  (3)  to  the  assigned  dominant  data  of  G,(s)  in 
(2),  we  can  directly  apply  the  dominant-data  matching  method  to  fit  a  low-order  stabilization 
filter  that  satisfies  the  specifications  assigned  in  (4).  Assume  the  desired  low-order  model  is 

From  the  definition  of  (4a)  we  have  b0  =  1.6a0.  and  (17a)  can  be  rewritten  as 

<17b> 

When  s  =  ju>„  =  /1.9  the  respective  values  of  | Ts2(J<tt„)l  and /t-,,  ;(/&*„)  in  (17b)  match  the  values 
of  |T,(ya>„)|  and/T,C/(o„)  in  (4c).  or  the  corresponding  non-linear  equations  are 

/,(a0,  a,.  b,)  =  2.56a0J  +  3.61h,2  -  2.576368889[(a0  -  3.61)2  +  3.61  n,2]  =  0  (18a) 


Ma„.  u„  b,)  =  1 .9  b,(a0  -  3.61)  -  3.04a0a,  -  0.075996381  ![I  .ba^ac  -  3.61) 

+  3.61a, b,]  =  0.  (18b) 

When  s  =  jwK  =  /'3.2  the  value  of/r,2(/wff)  in  (17b)  matches  the  value  of/r.Oow)  in  (4e),  or  the 
nonlinear  equation  is  ~  ““  “ 

/j(a o.  a„b,)  =  3.2b, (a#  -  10.24)  -  5. 12aoa,  -  0. 1 2798497821 1 .6 a<(a0  - 10.24) 

+  10.24a, b,]  =  0.  (18c) 

Applying  the  Newton-Raphson  method  [12]  and  using  the  initial  estimates  obtained  in  (15)  we 
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have  the  desired  solution  of  the  nonlinear  equations  of  (18)  at  the  9th  iteration  with  the  error 
tolerance  of  KT6.  The  desired  low-order  stabilization  filter  is 


T  856.6285965  +  21283. 19886  Q 

r2  +  3.318051s  +  13301.999297"  (  ' 

The  unit-step  response  curves  of  the  existing  stabilized  pitch  control  system  T,(s)  in  (la)  and 
the  redesigned  system  that  uses  the  low-order  filter  Tsi(s)  in  (19)  and  the  G0(s)  in  (le)  are 
shown  in  Fig.  3.  The  result  is  practically  identical.  From  the  response  curves  in  Fig.  3  we 
observe  that  T,2(s)  in  (19)  is  a  better  filter  than  the  T,,(s)  in  (16)  as  far  as  duplicating  the 
performance  of  the  original  pitch  control  system  is  concerned.  This  implies  that  the  existing 
stabilization  filter  Ts(s)  in  (lc)  might  be  overdesigned  since  we  can  duplicate  performance  with 
a  lower  filter.  Obviously,  the  implementation  cost  of  the  filter  Ts2(s)  is  less  than  that  of  Ts(s)  in 
(lc). 


HI.  AN  ALGEBRAIC  METHOD 

The  original  fourth-order  stabilization  filter  T,(s)  may  be  replaced  by  two  second-order 
filters,  T,i(s)  and  Ts,(s),  using  the  dominant-data  matching  method.  It  is  observed  that  all  three 
stabilization  filters  have  complex  roots  and  that  all  are  placed  in  the  feed  forward  loop. 
Therefore,  they  are  sensitive  to  external  disturbances.  If  alternate  filters  can  be  designed  and 
placed  in  both  the  feed  forward  and  feedback  loops,  then  (a)  the  designed  filters  may  result  in 
simple  transfer  functions  with  positive  real  roots  that  may  be  easily  synthesized  using  RC  type 
passive  elements,  and  (b)  the  reliability  and  cost  of  the  designed  system  may  be  improved.  The 
compensators  in  the  feedback  loop  enable  the  designed  system  to  be  less  sensitive  to  parameter 
variations  and  modeling  errors.  In  addition,  it  will  reduce  the  effects  of  many  external 
disturbances  [14]. 

The  algebraic  method  given  by  Shieh  [3]  and  Chen  [4]  is  extended  and  modified  to  redesign 
this  pitch  control  system.  The  steps  of  the  algebraic  method  are  summarized  as  follows: 

Step  1 

Assign  the  design  goals  using  frequency-domain  specifications  and  model  a  standard  transfer 
function  using  the  dominant-data  matching  method. 

Step  2 

Expand  the  obtained  standard  transfer  function  into  the  continued  fraction  expansion  shown 
in  (7)  to  obtain  the  dominant  quotients  and  to  formulate  the  matrix  equation  in  (9). 

Step  3 

Assume  the  fixed  configuration  compensators  with  unknown  parameters  and  determine  the 
overall  transfer  function  that  consists  of  the  unknown  parameters. 

Step  4 

Substitute  the  coefficients  of  the  obtained  over-all  transfer  function  in  Step  3  into  the  vectors 
[«]  and  [ft]  in  (9)  and  expand  the  matrix  equation  to  obtain  a  set  of  equations. 

Step  5 

Solve  the  set  of  equations  to  determine  the  unknown  constants  assigned  in  the  compen¬ 
sators. 

The  designed  system  using  the  algebraic  method  has  the  exact  dominant  quotients  of  the 
standard  models  and  is  a  good  approximation  of  the  standard  model. 

Before  we  design  the  pitch  control  system  using  the  algebraic  method  we  apply  the 
dominant-data  matching  method  to  determine  a  reduced-order  model  of  the  original  unstable 
system  G^s)  in  (te)  which  had  a  transfer  function  with  large  coefficients.  The  unstable  transfer 
function  Ggfr)  in  (le)  can  be  decomposed  into  a  stable  portion  and  an  unstable  portion  as 
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follows 


Go(j)  = 


I 


s(s~  2.921) 


Us) 


(20a) 


where  the  stable  portion  is 

T,  .  324332.316(s  +  0.1933X5  +  65X5  +  1500) 

*s)  (s  +  3.175Ks  +  87.9±/95.5Xi  +  1 12.5K-V  +  1385)' 

The  pole  at  the  origin  and  the  unstable  pole  at  s  =  2.921  are  considered  as  dominant  poles. 
Therefore,  they  are  retained  in  the  simplified  model  G0*(s),  or 

GM  =  Gf(s)  =  ,77^  290  ™  (20c) 

where  TS(s)  is  the  reduced  order  model  of  Us)  using  the  dominant-data  matching  method.  The 
frequency-response  data  of  the  gain  margin,  phase  margin,  phase-crossover  frequency,  gain- 
crossover  frequency  and  the  final  value  at  <a  =  0  are  used  as  the  dominant  data  for  the  transfer 
function  fitting.  The  resulting  T%(s)  is 

Tf.  .  496.854897s2  +192897.96101 1  s  +  37103.33375 

5 3  +  1 17.073733 +  16552.300003s  +  50595.685093'  (  ' 

The  desired  Tf(s)  is  a  low-order  model  with  small  coefficients.  Thus,  the  design  process  has 
been  greatly  simplified. 

Following  the  steps  in  the  algebraic  method  we  assign  the  series  compensator  C,(s)  and  the 
parallel  compensator  G^s)  with  unknown  parameters  Xhi  =  1, 2, ...  ,7  as 


C,(s)  = 


X,s  +  X 7 
s  +  X5 


(21a) 


and 


G2(s)  = 


X2s*+XAs  +  X2 
s2  +  X,s  +  X2  ' 


(21b) 


The  block  diagram  of  this  redesigned  system  is  shown  in  Fig.  4(1).  The  overall  transfer  function 


Fig.  4.  The  block  diagrams  of  (he  redesigned  system  using  algebraic  method. 


Computer-aided  methods  for  redesigning  the  stabilized  pitch  control  system 
Tf{s)  of  this  feedback  system  is 

T ,  v  ftp  +  bys  +  •  •  ■  +  A7S* 

'  S  <jp  +  a1s  +  ■  ■  •  +  ais* 
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where 


a0  =  37 103.33375  X2X7 

a,  =  192897.96101 1  X2X7  +  37I03.33375(X2X»  +  X4X7)  -  147789.996  »A';A's 
a,  =  496.854897X2X7  +  192897.96101 1(X2X«,  +  X4X7)  +  37103.33375(X4X* 

+  X,X7)  +  2246.4 1 679XjXs -  I47789.9961(X,  +  X,X,) 
a,  =  496.854897(X2X*  +  X4X7)+  192897.96101 1(X4X,,+  X,X7)  + 37103.33375 
X,X6-  147789.9961(X,  +  X,)  +  2246.41679(X:  +  X,X,)+  162 10. 32763  X2X, 
a4  =  496.854897(X„X<1  +  X,X7)  +  192897.96101 1X.X*  -  147789.9961  +  2246.41673 
•  (X,  +  X,)  +  16210.32763(X2  +  X,X,)  +  1 14.152733X,X, 

а,  =  496.854897X,X«,+  2246.41679  +  I62I0.32763(X,  +  X,) 

+  I14.I52733(X2+X,X,)+X2X, 
a„  =  16210.32763  +  1 14.152733(X,  +  X,)  +  X2  +  X,X, 
a7=  1 14.152733 +  X,  +  X, 
a»=  I 

b0  =  37 103.33375 X2X7 

б,  =  192897.96101 1X2X7  +  37103.33375(X;X<>  +  X,X7) 

b:  =  4%.854897X2X7  +  192897.96101  l(X2X*  +  X,X7)  +  37I03.33375(X,X*  +  X7) 
A,  =  496.854897(X2X*  +  X,X7)  +  192897.96101  l(X,Xk  +  X7)  +  37 103.33375  X„ 
bt  =  496.854897(X|X»  +  X7)  +  I92897.96I0I  I  X„ 

6,  =  496.854897 X* 

A«.=  0 

bj  =  0. 


In  order  to  match  the  seven  unknown  parameters  X,  *n  (21)  for  this  type  “1"  system  we  need 
eight  quotients,  hi  in  (9).  Therefore,  the  third  order  standard  model  in  (8b)  that  has  six 
quotients:  h,  =  I.  A2  =  -0.631845015,  A,  =  -0.476189214.  h,  =  14.799589050.  A,  =  -0.102867450. 
A*  =  -  13.924278040  should  be  increased  to  a  fourth  order  model  by  inserting  A7  =  100  and 
A*  =  0.1  where  the  values  of  A7  and  h,  are  selected  based  on  rules  developed  in(!5|.  The 
resulting  amplified  model  is 


Ti  \  637807  j 2055661s  +  0.243466s * 

As} "  6.37807  +  10.46222s  +  1.259008s3  +  s’ 

=  I 


A,  +  - 


A.  +  • 


A;  +  - 


A>  +  - 


A,  +  - 


A4  + 


A,  +  ~ 


^5  +’ 


A*  +  - 


^7  + 


K 


-.  63.78098007  +  21 1.8989926s  +  22.87561717s1  +  0.34346s* 

“  ' “ 63.78098007  +  1 10.9545225*  +  23.00917551^  +  1 1.301 105I5*3  +  s4'  (22) 


It  (1 5)  has  been  shown  that  (22)  is  a  good  approximation  of  the  original  model  in  (8b). 

Substituting  the  ah  i  =  0, 1 . 7  and  b,.i  =  0 . 7  in  (21)  and  using  the  quotients  hh  i  *  1 . 8 

in  (22)  into  (9)  yields  the  following  set  of  high-order  nonlinear  simultaneous  algebraic  equations. 

/,(X . X7)  =  X2X7  +  0.63 1 8422396[  X,(Xt  -  X,)  -  3.98319992X2X,1  -  0  (23a) 


I 

t 


m 
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/;(X . X7)  =  X^S. 22822291 X, + 8.522553 1 36X4  -  6.939879587X,  +  X,  -  I ) 

+  X}(I.582676549X„  -  13.17807554X0  +  X6(X4  -  X,)  t 

-3.983I99922(X:  +  X,X,)  =  0  (23b) 

/,(X . X7)  =  X2(  -  I2.7I36162IX*-  X,)+  X,(  13.58355291  X,  +  I.820964317X, 

-  26.29716913X0  +  !0.79844539(X1Xh  +  X7)  -  l3.31248704(X4Xk 
+  X,X7)  +  6.327224282(X,  +  X,)  +  1 .588477708X„<  I  -  X,) 

+  20.03527 1 43(X:  +  X,X,)  =  0  (23c) 

/4(X . X7)  =  X,(Xj  +  668.467007 1 X4  -  28 1 .48094X,)  +  X„(386.9860673Xj 

+  362.767005  -  456.258273X,)  -  647.2403649(X4Xfc  +  X,X7) 

-  57.53603068XjX,  -  548.5 188427(X2  +  X,X0  +  235.861385(X,Xk 

+  X7)  +  235.2945 185  +  590.5096275(X,  +  X,)  =  0  (  23d) 

/,(X . X7)  =  2357.408023(X,X*+  X7)  +  X„(  1 598.83993  IX,  +  17096.15228 

-  32881.95043X0  +  X,(4. 16745091 X,  +  1599.83991X,  -  X4) 

-  472.6735322(X4X*  +  X,X7)  +  24996.98242  -  939.0287936(XI 

+  X,)  -  2765.323026(Xj  +  X,X,)  -  52.0777 1943X,X,  =  0  (23e) 

/OX, . X7)  =  X„(  -  99.42094 1 5 X,  +  1 1 132.91981  X,  -  57256.87822)  +  XAX* 

-  1 00.420941 5 X,)  +  41 1 .4274907(X4X„  +  X,X7)  +  67006.93001(X, 

-  xo  +  1234.567433(Xj  +  X,XO  +  42.6901 1 171X2X,  +  23203.53455 

-  391 12.69694(X,X*  +  X7)  =  0  (23f) 

MX . X7)  =  4%.854897(XiX*  +  X,X7)  +  198512.9704(XIX»  +  X7) 

+  2228495 .695  X*  -  1 70.649783 1(X4X*  +  X,X7)  -  77618.59617X,Xk 

-  3442861.087  -  395845.4335(X,  +  X,)  -  8390.81 2346(X«  +  X.XO 

-  62.0825 1489X:X,  =  0.  (23g) 


Equation  (23)  is  a  high  order  nonlinear  equation  which  is  difficult  to  solve.  The  Newton- 
Raphson  method!  12)  is  applied  to  determine  the  unknown  parameters  X,.  The  following  method 
is  suggested  for  obtaining  the  initial  estimates. 

The  structure  of  the  desired  fixed  configuration  control  system  in  Fig.  4(1)  can  be  modified 
as  shown  in  Fig.  4(2).  The  overall  transfer  function  is 


(24a) 


where 


.  Gi(s)G2(s)Gj(s) 

,jW“|  +  G,(s>Gj(s)GS(s) 

The  design  objective  is  to  determine  G,(s)  and  Gj(s)  such  that  the  response  of  r,(5)  is  close  to 
that  of  the  standard  model  T,(r)  in  (8b).  Replacing  the  series  compensator  G,(s)Gj(s)  in  Fig. 
4(2)  by  the  designed  stabilization  filter  T,£s)  in  (19)  and  equating  the  resulting  transfer  function 
T,(s)  in  (24a)  to  the  standard  mode!  T,(s)  in  (8b),  we  can  solve  the  approximate  transfer 
function  Gt(s)  of  G7(j),  or 


G?(j)  =  G^s) 


IM 

T,(i) 


GMGMGtU) 

1  +  G,(j)G^i)GS(j)rf(j) 


TM*)Gt(s) 

I  +  TMs)Gtb)Trb) 


-  [s.036619205  x  I0»+  3.46495752  x  I0'°s  +  4.540060393  x  10'V  +  7.840679235  x  10’s' 

+  4.363076841  x  l0*s*  +  1.763524302  x  I0V  +  4.256201 128  x  I05j‘)/[5.0366I9205 
x  10*  +  3.008227329  x  I0,#s +  4.613716606  x  IO'V  +  6.124169121  x  10V 
+  4.498497844  x  10’s4  +  8.512494768  x  10V  +  9.985459768  x  lO's*  +  9.698650697  x  10’s7 
+  4.915681 19  x  I0s'  + 0.243466s*].  (24b) 
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Expanding  (24b)  into  the  form  of  (7a)  yields  a  set  of  dominant  quotients  A,  =  1,  A2  = 
-  1.102755917,  A,  =  -0.1287948973,  ht  =  5.593229805.  A,  =  0.1338916858.  A«  = 

•  •  •.  Substituting  the  first  five  quotients  into  (7e)  gives  a  second  order  approximate  model 
G2**(s)  of  the  approximate  parallel  filter  G$(s)  in  (24b)  yielding 


G?*(f)  = 


0.994929057s2  +  0.7394973923s  -I-  0.1058245527 
s1  +  0.643533679s  +  0.1058245527 


(24c) 


Ct*(j)  is  an  approximate  model  of  the  assigned  parallel  compensator  Gj(s)  in  (21b).  The 
approximate  series  compensator  (Gt(s))  of  G,(s)  is 

T,M  2252.284999+  13787.1076s  +  21834.46821s2 
T(s'  G!*(s)  1407.678125  +  9837.144919s  +  13237.10512s2 

+  856.628596s* 

+  4.040722745  s*  +  0.994929057s*' 


Equation  (25a)  can  be  expanded  into  the  form  of  (7a)  to  obtain  a  set  of  dominant  quotients 
A,  =  0.625,  A2  =  1.845828612.  A,  =  0.0839039052,  A4  =•••,•••,  A,  =  ••  •.  Substituting  the  first 
three  quotients  into  (7c)  yields  the  reduced-order  model  of  Gf(s)  in  (25a) 


Gt*(s)  = 


1.410628426s +0.2184671685 
s  +  0.1 3654 19803 


(25b) 


Gt*(s)  is  an  approximate  model  of  the  assigned  series  compensator  G,(s)  in  (21a).  Comparing 
(21b)  and  (24c)  and  (21a)  and  (25b)  we  have  a  set  of  initial  estimates  as  Xf  =  0.643533679,  XJ  = 
0.1058245527,  X\  =0.994929057.  XJ  =  0.7394973923,  X?  =0.1365419803,  XJ  =  1.410628426  and 
X?  =  0.2184671685.  Using  these  initial  estimates  and  the  Newton-Raphson  method[12]  to  solve 
the  nonlinear  simultaneous  algebraic  equations  in  (23)  yields  the  solution  X,  =0.503850, 
X2  =  0.059928,  X,  =  1.051503,  Xt  =  0.580016,  X,  =  4.831826,  X6  =  1.885577,  and  X7  =  6.744450  at 
the  14th  iteration  with  the  error  tolerance  of  10”*.  The  desired  compensators  are 


G,(s)  = 


1.885577s  +6.744450  1.88SS77(s  +  3.57688) 
s +4.831826  s +4.831826 


(26a) 


and 


l.05!503sJ  +  0.5800I6s  +  0.059928 
s2  + 0.503850s +  0.059928 
1.051503(s  +0.13769KS  +0.41391) 
(s+O.I9244Ks+0.3U405) 


(26b) 


The  unit-step  response  curves  of  the  existing  stabilized  system  T,(s)  in  (la)  and  the  redesigned 
system  using  the  compensators  in  (26)  and  the  G«(s)  in  (te)  are  shown  in  Pig.  3.  The  result  is 
satisfactory.  Note  that  G,(s)  and  Gj(s)  in  (26)  are  positive  real  functions  with  positive  real  poles 
and  zeros  that  can  be  realized  using  RC  type  passive  elements.  Thus,  the  advantage  of  a 
feedback  controller  system  [14]  may  be  fully  utilized. 

IV.  CONCLUS’ON 

Two  computer-oriented  methods:  a  dominant-data  matching  method  and  an  algebraic 
method  have  been  presented  to  redesign  an  existing  stabilized  pitch  control  system  [  I  ].  Thus,  an 
alternate  method  to  the  trial-and-error  approach  that  is  traditionally  used  has  been  given.  The 
resulting  low-order  stabilization  fitters  that  were  obtained  using  the  above  methods  reduce  the 
implementation  cost  of  the  pitch  control  of  the  missile  system. 

The  dominant-data  matching  method  can  be  applied  to  general  control  system  design  to 
determine  the  reduced-order  model  of  a  high-order  system.  The  proposed  model  reduction 
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method  is  superior  to  most  existing  methods  in  that  it  provides  the  exact  frequency-domain 
specifications.  I 

An  algebraic  method  has  been  applied  to  determine  fixed  configuration  filters  such  that  the 
performances  of  the  redesigned  pitch  control  system  of  an  example  missile  has  been  greatly 
improved.  The  algebraic  method  can  be  applied  to  design  any  desired  control  system  structure. 

Several  methods  have  been  given  for  estimating  good  initial  guesses  for  solving  nonlinear 
equations.  In  summary,  the  proposed  computer-aided  design  methods  present  an  attractive 
alternative  to  trial-and-error  methods  which  can  be  used  to  design  control  systems. 
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